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Abstract 

The HYPERDIRE is a project devoted to creation of a set of Mathematica based 
programs for differential reduction of hypergeometric functions. The current version 
inchides two parts: one, pfq, is relevant to manipulation with hypergeometric func- 
tions p4_iFp, and second, AppellFlF4, to manipulations with Appell's hypergeometric 
functions of two variables Fi, F2, F^, F4. As application, a few multiloop Feynman Di- 
agrams of propagator type are considered. 
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LONG WRITE-UP 



1 Introduction 

Multiple hypergeometric functions [lrt4j play an important role in many branches of sci- 
ence. In particular, a large class of Feynman Diagrams are expressed in terms of Horn-type 
hypergeometric functions [5]. 

Let us consider a multiple series: 



H{f, <?; f ) = 



with flab, i^ab G Q, 7i, CTfc G C. The sequences 7 = (71, ■ ■ ■ , 7^-) and a = (ai, ■ ■ ■ , 0"^) are 
called upper and lower parameters of the hypergeometric function, respectively. 

Let Cj = (0, ■ ■ ■ , 0, 1, 0, ■ ■ • , 0) denote the unit vector with unity in its j^^ entry, and 
define = ■ ■ ■x^'' for any integer multi-index m = (mi, ■ ■ • ,771^). Two functions of 
type (II]) with sets of parameters shifted by unity, + e^; a; x) and H(j] a; of), are related 
by a linear differential operator: 



+ ec;a;x) = — \Y1 f^caXa^ + 7c ) H{f, a; x) = f/[+ ^^^^yiJ(7, a, x) 

\a=l "^'^ J 



(2) 



Similar relations also exist for the lower parameters: 

H{l\ o - e^- x) = — ^ + j ^' ^) = ^[a,^a,-i]H{f, o"; x) . (3) 

The linear differential operators ?7^_>.^^_,_i, L~^^^^_-^ are called the step-up and step-down 
operators for the upper and lower index, respectively. If additional step-down and step-up 
operators f/~, satisfying 

f/r , 1^ ^U^ ^ , (7, a, x) = L,^ ]L7 ^ ^.H(^,a,x) = HH,cr,x) 

[7c + l— S-7c] |7c— ^7c + lJ ^ " ' [(Jc — 1— S-(Tc] ^CTc— >CTc-lJ \ I 1 1 J \ I 1 7 J 

{i.e., the inverses of U^^, L~J are constructed, we can combine these operators to shift the 
parameters of the hypergeometric function by any integer. This process of applying t/^, L^^ 
to shift the parameters by an integer is called differential reduction of a hypergeometric 
function. 

In this way, the Horn-type structure provides an opportunity to reduce hypergeometric 
functions to a set of basis functions with parameters differing from the original values by 
integer shifts: 

E I'll / (9 \ ™^ / (9 \ 

Po{x)H{^+k;a + l-x)= Yl Pmu-,mM 1^-] ■■■[g-j H{f,ff;x), (4) 



where Poix) and Pmi,---,mp(a^) are polynomials with respect to 7, a and x, and k, I are lists of 
integers. 

Algebraic relations between the functions -^(7, <?; x) with parameters shifted by integers 
are called contiguous relations. The development of systematic techniques for solution of 
contiguous relations has a long story. It was started by Gauss who described the reduction 
for the 2F1 hypergeometric function in 1823 Numerous papers have since been published 
[OHH] on this problem. An algorithmic solution was found by Takayama in Ref. [9j, and those 
methods have been extended in a later series of publications [101111] (see also |T2]). 

Let us recall that any hypergeometric function can be considered to be the solution of 
a proper system of partial differential equations (PDEs). In particular, for a Horn- type hy- 
pergeometric function, the system of PDEs can be derived from the coefficients of the series. 
In this case, the ratio of two coefficients can be represented as a ratio of two polynomials. 



C{m + ej) _ P.j{m) 



C{m) Qji'rn) 



(5) 



so that the Horn-type hypergeometric function satisfies the following system of differential 
equations: 



= A' (7, ^,x)H{j, a,x) 




H{^,a,x), (6) 



where j = 1, ■ ■ ■ , r. It was pointed out in several publications [T3HT5] that (i) the differential 
reduction algorithm, Eq. @, can be applied to reduction of Feynman Diagrams to some 
subset of basis hypergeometric functions with well-known analytical properties p^[T^ : (ii) 
the system of differential equations, Eq. (E]), can be used also for construction of, so-called, 
£-expansion of hypergeometric functions around rational values of parameters via direct 
solution of the system of differential equations |T5] . This is another motivation for creating 
the package for manipulation with parameters of Horn-type hypergeometric functions. 

The aim of this paper is to present the Mathematica [16] based package HYPERDIRE 
for the differential reduction of the Horn-type hypergeometric function with arbitrary values 
of parameters to a set of basis functions. The current version consists of two parts: one, pfq, 
for manipulation with hypergeometric function, p+i-Fp, and the second one, AppellFlF4, for 
manipulations with Appell functions, Fi, F2, F3, F4. The algorithm of differential reduction 
for other functions can be implemented as an additive module. 

In contrast to the recent programs written by members of Computational Particles 
Physics [T7I - [TU] . the aim of our package is the manipulation with hypergeometric functions 
without construction of ^-expansion [20l|2T] . 

The preliminary version of pfq was presented in [22] and is available in [23]. The latest 
version is available in 
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2 Differential reduction algorithm for generalized hy- 
pergeometric function p+iFp 

2.1 General consideration 

Let us consider the generalized hypergeometric function, pFq{a; b; z) , defined around z = 
by a series 

,f,(„-s;.).,n(^:|.)^g£;5y^|. (7) 

where {a)k is a Pochhammer symbol, {a)k = r(a + k)/r{a). The sequences a = (ai, ■ ■ ■ , Op) 
and b = {bi, ■ ■ ■ ,bq) are called the upper and lower parameters of hypergeometric functions, 
respectively. In terms of the operator 6: 

the differential equation for the hypergeometric function pFq can be written as 

[zUtiie + a,) - euUiO + h- l)] pF,(a; b; z) = 0. (9) 

2.2 Differential reduction 

The Differential reduction for these functions was analyzed in details in |T^. Here we remind 
some of the main relations relevant to our program. 

The universal differential operators, Eqs. ([2]), ([3]), have the following form: 

pFq{ai + l,d]b]z) = B^^pFg{ai,d;b;z) = — {9 + ai)pFg{ai,d;b;z), (10) 

pFq{d;bi-l,b;z) = H.pFg{d;bi,b; z) = [O+bi-l) pFq{d;bi,b] z) , (11) 

Oi — i 

where the operators B^^{H^) are called the step- up (step-down) operators for the upper 
(lower) parameters of hypergeometric functions. This type of operators were explicitly con- 
structed for the hypergeometric function p+i-Fp by Takayama in Ref. [lOj. For completeness, 
we reproduce his result here: 



— * — * 

.iFp{ai - 1, a; b; z) = B~ p+iFp{ai, a; b; z) , 

yiFp{d]bi + l,b]z) = H^, p+iFp{d]bi,b; z) , (12) 
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where 



b: 



Uix) = 



xlVj^i{x + bj -1) - Ci 



X + Qi 



\x- 



i+1 . 



j=0 



X + Qi 



j=0 



(13) 



fe=0 



and 



H. 



+ 



bj-l 
di 



di = n^;i'(i + a, -6,) , 

Si{x) = 



Ii%\{x + a,)-d: 



p(P+i) / r„ Ax'-jl-hiY] 



3=0 



(14) 



fc=0 



where 
as 



(f) is a polynomial of an order of s with respect to variables r, and 

iV=l ii<---<ij 



means substitution of a by a + 1, and polynomials Pj^\ri, ■ ■ ■ ,rp) are defined 



P 



i=o j=0 j=Q 



(15) 



fe=l 



For example, P'f\r) — Y7j=i '^j P^\r) — 11^=1 '^j- Keeping in mind that 

p p+k p+k 

we get that these polynomials satisfy the following relations: 

k 

P^%%{n, ■■■ ,rp,qi,--- ,qk) = Yl ^P+i-.-n(n, • • • , r,)Pi'\q^, ■■■,q,), (16) 

n=0 



6 



where j = 1, ■ ■ ■ ,p — k. In particular, 

P',l\'},{r,_i.qi.q2) = Y.Pt\'},-u{r)P^k\<i) ^ j = l,---,p-l. 



fc=0 
3 



k=0 



The differential reduction has the form of a product of several differential step-up/step- 



down operators, and , respectively: 



Fia + m;b + n;z) = ( Hf^^ 



Bf ^ 



{b} ) ^(«; b; z) , 



(17) 



so that maximal power of 9 in this expression is equal to r = Yli'^^i + Xlj'^i- Since the 
hypergeometric function p+iFp(a; b; z) satisfies the differential equation of an order of p + 1 
(see Eq. (ED): 

{l-z)e^^\+^F^{a-h-z) 

EkitUK})-^i?i-.({&.-i}) 



e' + zUX:Xak)p+iFp{d-h-z), (18) 



, r=l 



it is possible to express all terms containing higher power of the operator 6'^, where k > p + 1, 
in terms of rational functions of parameters and argument z times by 9\ where j < p. 

In this way, any function p_(.iFp(a + m; h+k\ z) is expressible in terms of the basic function 
and its first p-derivative: 



p+iFp(a + m; h + k]z 
1 



(19) 



S{ai,hj,z) 



Ri{ai, bj, z) + Riidi, bj, z)6 -\ h Rp+i{ai, bj, z)9^ L+iFp(a; 6; z) 



where m, k is the set of integer numbers, and S and Ri are polynomials in the parameters 
{oi}, {bj} and z. 

From Eq. f|T3|) follows that if one of the upper parameters aj is equal to unity, then 
the application of the step-down operator B^^ to the hypergeometric function p+i-Fp will 

produce unity, p+iFp(l, a; &; 2) = 1 . Taking into account the explicit form of the step- 
down operator B^ , 



1 



K=iih - 1) 



(20) 
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we get the differential identity 

Yi%, {bj - 1) -zn^^ia, + {l~z)9P^p+^Fp{l, a; b; z) 
■p-i >, 

E - 1}) e^yiFAi,d-b-z) = nLi(&.-i) . (21) 

The case when two or more upper parameters are equal to unity, Oi = 02 = 1, does not 
generate any new identities. 

3 Differential reduction for the Appell hypergeometric 
functions 

3.1 Appell hypergeometric functions: system of differential equa- 
tions 

Let us consider the system of linear differential equations of the second order for the functions 

uj{z): 

eiM^ = |Po(^)^l2 + Pl(^)^l+P2(^1^2 + P3(^)|c^(^1 , (22) 
^22W(£) = |i?o(?)^12 + i?l(^)^^l + i?2(^)^2 + i?3(^)|c^(^) , (23) 

where z = (^1,^2) and zi,Z2 are variables, {Pj,i?j} are rational functions, 9j = zjdz for 
j = 1,2 and ^ji-.-i^. = 6i.- ■ -Oi^. The exploring of 9j instead of the standard dj is explained 
by our applications. Taking the derivative of Eq. fj22|) with respect to ^2? using the well- 
known property d2diiuj{z) = didi2Uj{z) and applying Eq. fl23|) . we rewrite the first equation 
as follows: 

[Ol-Po02] 012i^{z) 

[92P0 + P1 + P2M ei2+[P2Rl + e2Pl] 9i + [P2R2+92P2 + P3\ 02 + P2R3 + 02P3^CO {z} , 

(24) 

Repeating a similar operation for Eq. ( 123|) we get 

[-Roei+e2] ei2uj{z) 

[eiRo+R2+RiPo] 9i2 + [PiRi+0iRi+R3]ei + [P2Ri+9iR2]92+RiP3+9iR3^uj{z) , 

(25) 
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It is well known [2] that under condition 

1 - PoRo ^ , (26) 

there are 4 independent solutions of system fl221) . fl2^ . In this case, Eqs. . fl2S]) can be 
solved so that 
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{l-PoRo)en2Co{^ = YPoi0iRo + R2 + RiPo) + 02Po + Pi + P2Ro]O 

+ [P^R, + 92P1 + Po {PiRi + e^R^ + Rs)] 9^ 

+ [P2R2 + e2P2 + P3 + Po (P2i?l+^li?2)] ^2 

+P2R3+e2P3+Po {RiPs+e.Rs) \uj{z) , (27) 



{l-PoRo)9i22Co{^ = Y^iRo+R2+RiPo+Ro{02Po+Pi+P2Ro)]0l2 

+[Ro {P2Ri+e2Pi)+PiRi+e,R^+R3\ Oi 
+[P2Ri+eiR2+Ro {P2R2+O2P2+P3)] O2 

+Ro {P2R3+92P3) + RiP3+9iRs^u{z) . (28) 

The condition of complete integrability is defined from the relation di {di22^^{z)) = 82 (9ii2w(i')), 
or in terms of the quantities ( 1271) and ( l28l) it has the following form: 

= )rS- [r-h.s.Eq. m] - -, ^ |r.h.s.Eq. gai 

+e2{e^2^{z))-e^{e^2^{z)) , (29) 

where we have used the following relation 9jj = z|(9| + Zjdj. If the condition f l2I?]) is valid, 
Eqs. (122|) and (123|1 can be reduced to the Pfaff system of four differential equations 

df = Rf, (30) 

where / = {ijj{z),9iOj{z),92ijj{z),9l2'-^{z)) . 
In the case 

1 - PoRo = , (31) 
9i2Uj{z) are expressible in terms of three other elements: u{z),9iu{z),92Uj{z). In particular. 
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Table 1: 



Rq 



Fi 



1 



(a+bi)zi~ 



7^ 



1-zi 

Z2 
1-22 



(a+bi)zi-(ci-l) 
^^£1^ 



(1-^1 ) 

" 1-22) 



(ai+fei)2i-(^^^ 
^^£1^ 



^4 

221 

(l-zl-22) 

222 

(1-21-22) 

(a+6)2l-(Cl-l)(l-22) 
(1-21 -22) 



Pi 



(1-21 



(g+b+l- £1)2:2 
(1-21 -22) 



Ri 



^2 22 



(1-22 



^222 







(g+b+l- C2)2i 
(1-21 -22) 
(a+6)22-(c2-l)(l-2i) 

(1-21 -22) 



P, 



biZi 



(1-21 



(0+62)22- 



bizi 
(I-21) 
(0+62)2:2 -(C2-l) 







(a2+fc2)22-(c-l) 





R2 



(1-22 



ab\z\ 



ab\z\ 
ab2Z2 



a\b\z\ " 



ab2i 



(1-21 

ab2 22 



a2fa2 22 
(1-^2) 



(1-21 -22) 

ab22 



(1-22 



(1-2:1-22) 



in terms of the notation of Eqs. 

[^2^0 + Pi + P2R0 + Po (O.Ro + R2 + RiPo)] Ouco{^ 

- [Po {OiRi + R3 + RiPi) + P2R1 + O2P1] Oi 

- [Po (^li?2 + R1P2) + P2R2 + ^2i'2 + P3] 02 

- [Po (e.Rs + P1P3) + P2R3 + 02P3] \tuiz) . 



(32) 



In this case, the integrabihty conditions are vahd and Eqs. (1221) and (!23|l can be reduced 

— * 

to the Pfaff system ( l30l) of three differential equations / = {u{z),6iu{z),92Uj{z)) , and the 
system has 3 solutions. 

For Appell's hypergeometric functions Fi,F2,F3 and F4 the values of coefficients of 
Eqs. ([22D, ([23D are collected in Tabled 



3.2 Appell hypergeometric function: Fi 
3.2.1 General consideration 

Let us consider the Appell hypergeometric functions Fi defined around x = y = as 

{a)m+n{bl)m{b2)nX"' 



u = Fi{a,bi,b2,c;x,y) = 



m=0 n=0 



)m+n 



(33) 



In this case, Eqs. (1221) and (1231) have the following form: 



0xx^ 



9yyU 



{a + bi)x-{c-l) y 
02 



1 — a; 



x-y 



^'^^^ - e..u + -^ahuj , (34) 



1 — a; 



{a + h2)y-{c- 1) , X 
h 61- 



x~y 



^^y^^-'^ 9^u + -^ab2U. (35) 
1 
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The Eq. ( 1311) is fulfilled and Eq. ( 1321) has the following form: 

, , d^uj ^ duo ^ duj 
oxoy ox oy 



or in terms of the operators Ox-, Oy 



OxyUJ 



hix 



■6yU 



X — y X — y 

3.2.2 Differential Reduction of Fi 

The direct differential expressions follow from Eqs. ([2]) and ([3]) 

aFi{a + l,bi,b2,c;x,y) = {9x + 9y + a)Fi{a,bi,b2,c;x,y) , 
biFi{a,bi + l,b2,c;x,y) = (6'^. + 6i)Fi(a, 6i, 62, c; x, , 
(c- l)Fi(a, 61,62, c- = {9cc + 9y + c-l)Fi{a,bi,b2,c;x,y) 

The inverse differential relations were considered in [2l[7]: 

(c-a)Fi(a - 1, 61, 62, c; x, y) = 

[c-a-bix-b2y+{l-x)9x + {l-y)9y] Fi(a, 61, 62, c; x, ?/) , 
(c-6i-62)Fi(a, 61 - l,b2,c;x,y) = 

c-bi-b2-ax+{l-x)9x-x ^1-^ 9y Fi{a,bi,b2,c;x,y) , 



(c-a)(c-6i-62)Fi(a, 61, 62, c + 1; x, y) 



(c-a-61-62 



9. 



9,, 



Fi{a,bi,b2,c;x,y) 



(36) 



(37) 



(38) 
(39) 
(40) 



The differential reduction for the parameter "62" follows from Eqs. (139 p and (142 p 
symmetry property of the function Fi. Fi{a, 61, 62, c; x, y) = Fi{a, 62, 61, c; y, x) , so 

bi^b2 , X . 



(41) 
(42) 

(43) 

and the 
that 



3.3 Appell hypergeometric function: F2 
3.3.1 General consideration 

Let us consider the Appell hypergeometric functions F2 defined around x = ?/ = as 

(a)m+n(&l)m(&2)nX"' ?/" 



uj = F2(a, 61, 62, ci, C2; x,y) = ^Y1 



(ci)m(c2)„ n! 



m=0 n=0 

In this case Eqs. ( I22l) and ( l23l) have the following form: 

(1 — x)^^2:2:Ci; = x9xyUJ+[{a+bi)x — {ci — l)]9xUJ+bix9yijj + abixu, 
{l-y)9yyU = y9xyUJ + [{a+b2)y-{c2-l)]9yUJ + b2y9xUJ + ab2yuj . 



(44) 



(45) 
(46) 
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The condition, Eq. (l26l) . is fulfilled and Eqs. ( 1271) and ( l28l) have the following form: 

h2xy 



(a + 6i + l-C2)a;-(ci-l)(l-i/) + 

I f^^^ [a + 6i + l-ci] 6^^;^ 
1 — a; 



1 — X 



+ 



(a + l-C2)&ix4 



1—x 



9yU 



abib2xy 

+ — u 

1 — x 



{l-x-y)9yyxU 



(a + 62 + l-Ci)i/-(c2-l)(l-x)4 



bixy 



1-1/. 



OxyUJ 



_^bixy^ [a+fe2 + l-C2] 
1-7/ 



bib2xy 



+ 



-1/ 

(a+l-ci)622/4 



abib2xy 
+ — UJ . 



3.3.2 Differential Reduction of F2 

The direct differential expressions follow from Eqs. ([2]) and ([3]) 

0^2(0+ 1,61, 62, ci,C2;x,?/) = {a+9cc + 9y)F2{a,bi,b2,ci,C2;x,y) 
biF2{a,bi + l,b2,ci,C2;x,y) = {bi + 9x)F2{a,bi,b2,ci,C2; x,y) , 
{ci-l)F2{a,bi,b2,Ci - l,C2;x,y) = (ci -l + 6'^.)^2(a, &2, Ci, C2; x, 

The inverse differential relations were considered in [7]: 

/ 1 L L ^ /1 xbi-{l-x)9x yb2-{l-y)9y 
F2{a - \,biMiC\,C2\x,y) = < 1 



Ci + C2-a-l 



1 1 
+ 

Ci — a C2 — a 



xF2(a,6i,62,ci,C2;x,?/) , 
^2(0, 61 - l,b2,ci,C2;x,y) 

F2{a,bi,b2,ci + l,C2;x,y) 



Ci — a C2 — a 

[( 1 - a; - y ) 6*^.5^ - 6ia;6'y - 62^/6' J 



(47) 



(4J 



(49) 
(50) 
(51) 



(52) 

'2,ci,C2;x,y) , (53) 



(ci-a)(ci-6i) "I ^ ( X ' 



1 



x(ci + C2 — a — 1) 



[xbi9y+yb29x-{l-x-y)9a:y] \ F2{a,bi,b2, ci, C2; x,y) . (54) 
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The differential reduction for tlie parameters "62" and "C2" follows from Eqs. (1501) . (1531) and 
Eqs. (l5Ti) . (!5^ and symmetry property of function F2: -^2(0, &2, c; y) = ^2(0, 62, c; |/, x) 

61 62 , ci C2 , X . 

3.4 Appell hypergeometric function: F3 

3.4.1 General consideration 

Let us consider the Appell hypergeometric functions F3 defined around x = y = as 

u = Fsia^, a2, h, 62, c;x,y) = }^}_^ — — , (55) 

m=0 n=0 

In this case Eqs. f l22p and (12 3 p have the following form: 

{l — x)9xxOJ = —9xyLo + [{ai + bi)x — {c—l)]9xU+xaibiU, (56) 
{l-y)6yyUJ = -6:ryU} + [{a2 + b2)y-{c-l)]6yUj + ya2b2UJ . (57) 

The condition, Eq. fl26|) . is fulfilled and Eqs. fl27|) andf l28|) have the following form: 

{xy-x-y)9^^yU = [{l-y){ai + bi)x-y{a2 + b2 + l-c)]9^yUJ 

+ {l — y)xaibiOyUj — ya2b26xUJ , (58) 

{xy-x-y)6^yyU = [{l-x){a2 + b2)y-x{ai + bi + l-c)]6^yUJ 

+ {l — x)ya2b20xOJ — xaibi6yUJ . (59) 

3.4.2 Differential Reduction of F3 

The direct differential expressions follow from Eqs. ([2]) and ([3]) 

01^3(01 + l,a2,&i,&2,c;x,?/) = (ai + 6'^)F3(ai,a2,&i,&2,c;x,y) , (60) 
61^3(01,02,61 + 1,62, c;x,?/) = (6i + 6'^.)-^3(ai, 02,61, 62, c; X, y) , (61) 
(c-l)F3(ai,a2,,6i,62,c- l;x,y) = (c-l + 6'^ + 6'y)F3(a, 61, 62, Ci, C2; x, y) . (62) 
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The inverse differential relations were considered in [8]: 
^3(01 - 1, 02, &i, &2, c; X, y) = 1 + 



{c-ai-a2){c-b2-ai) 



xF3(a, 61, 62, ci, C2; X, y) , (63) 

^3(01, 02, 61 - 1, &2, c; X, = 1 + — — 

[c-hi-h2)[c-a2-hi) 



yj \ y, 

xFsi^a, bi, 62, ci, C2; X, y) , (64) 
AF3(ai, 02, 61, 62, c + 1; X, ?/) = 

xF3(a, 61, 62, ci, C2; X, y) , (65) 
where in the last expression, 

61 = c—ai—bi , 62 = 0—02 — 62 ; F = 0—01 — 02 — 61 — 62 5 
A = (5i52-f'+«i6i5i +0262(^2 , 

Di = 52-^+0161-0262 , D2 = SiF+a2b2-aibi , B = Si + S2 , 

A = (c-6i-62)(c-oi-02)(c-02-6i)(c-oi-62) . (66) 

The differential reduction for the parameters "02" and "62" follows from Eqs. (!60l) . (163!) and 
Eqs. (16T|) . (IMI) . respectively, and the symmetry property of the function F3 : 
^3(01, 02, 61, 62, c; X, y) = ^3(02, Oi, 62, 61, c; y, x) : 

Oi <^ 02 , 61 <^ 62 , X 4^ y . 

3.5 Appell hypergeometric function: F4 
3.5.1 General consideration 

Let us consider the Appell hypergeometric functions F4 defined around x = y = as 



00 00 



The condition, Eq. ( !26l) . is fulfilled and Eqs. ( !27l) .(l28l) have the following form: 

{l — x — y)6xx^ = IxOxyOJ + [(o + 6)x — (ci — 1)(1— y)] OxOJ + [0 + 6+I — C2] x^^ya; + o6xa;(|38) 
{l — x — y)6yyUJ = 2y9xyUJ + [(0+6)?/— (c2 — 1)(1 — x)] + [a+b+l — Ci]y9xUj+abyijj (.69) 
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The condition, Eq. (l26l) . is fulfilled. However, instead of Eqs. (1271) . (l28l) . which are very long 
in this case, we present the results in the following form: 



[{l-x~yf-4:xy] O^^yU = 

y [l — x — y+2x{a+b+l — ci)] O^x^ + x [2x + (l — x — ?/)(a+6+l — C2)] OyyUj + 

2(a + 6+l-C2)(l-x)a; - {ci-l){l-x-yf - x{l-x-y){a + h+Ci-l) 
+y[{ci-l){l-x-y) + 2abx\exUJ + x[2x{c2-l) + {l-x-y)ab]eyUJ , (70) 



9,j.yUJ 



[{l~x-yf-Axy] 6. 



yyx 



X [l-x-y+2y{a+b+l-C2)] OyyU + y [2y+{l-x-y){a + b+l-ci)] O^x^ + 

2{a + b+l-ci){l-y)y - {c2-l){l-x-yf - y{l-x-y){a+b+C2-l) 
+x [(c2-l)(l-a;-y) + 2aby] 6yU + y [2y{ci-l) + {l-x-y)ab] 9xOJ , 
where the values of Oxx^^ and 9yyU are taken from Eqs. ( I68l) and ( l69l) . 



O^yUJ 



3.5.2 Differential Reduction of F4 

The direct differential expressions follow from Eqs. ([2]) and 



0^4(0 + l,b,ci,C2;x,y) 
(ci-l)F4(a, 6, ci - l,C2;x,y) 



{a+9^ + 9y)F^{a, b, ci, Ca; x, y) , 
(ci-l+6'a:.)-F4(a, 6, ci,C2;x,?/) , 



(71) 



(72) 
(73) 
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The inverse differential relations were considered in [7]: 

Fi{a-l,b,ci,C2;x,y) = |l 
1' 



X 



(ci-a) 



1-- ) 9.+9y+b 



1 



y 



(c2-a) 
1 



{1 — X — y){ci + C2 — a — 1) [ci — a C2 — a 



+ 



xy 



[(1 - a; - - Axy]e, 
-y[2x{a + b+l-ci) + {l-x-y){b+l-ci)]e^ 
-x[2y{a + b+l-C2) + {l-x-y){b+l-C2)]ey 



-2abxy 



>Fi{a,b, ci,C2;x,y) 



(74) 



x{ci-a){ci-b)F4{a,b,Ci + l,C2]x,y) = 
ci{ x{ci — a — b) + {l — x — y)6x — —{a + b+l — C2)9y 

Cl 



(2ci+C2-a-6-l) 



(l-x-t/)(ci + C2-a-l)(ci + C2-&-l) 



[{l-x- yf - Axy]e, 



xy 



-y [(l-x-?/)(a + 6-C2-2ci + 2) + 2x(a+6-ci + l)] 9^ 
-— [2ciy{a+b+l-C2) + (l-x-y)(c2-a-l)(c2-&-l)] Oy 

-2abxy ^'^F4^{a,b,ci,C2;x,y) . (75) 

The differential reduction for the parameters "6" and "C2" follows from Eqs. f l72|) . (17^ and 
Eqs. ( 173|) . fl75l) . respectively, and the symmetry property of the function F4 : 

(i) Fi{a,b,Ci,C2,c;x,y) = Fi{b,a,ci,C2,c;x,y) : a <^ 6 , 

(ii) F4(a,6,Ci,C2,c;x,?/) = F4(a,6,C2,Ci,c;?/,x) : ci <^ C2 , x^y. 
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Table 2: 





{a, 61, 62, c— a, c— 61 — 62} G ^ 


F2 


{a, 61, 62, ci-a,C2-a,ci + C2-a, 01-61,02-62} G ^ 


F3 


{ai, 02,61, 62, c-ai-a2, 0-61-62,0-02-61,0-01-62} G ^ 


F, 


{a,6,oi—a,oi — 6,02 — a, 02 — 6,01+02 — a, 01 + 02 — 6} G Z 



3.6 Appell hyper geometric functions: exceptional values of pa- 
rameters 

As was explained in Section ISTTl the differential reduction algorithm in application to Appell's 
functions could be written symbolically as follows: 

R{x,y)Fi{A + m;x,y) = [Po{x,y) + Pi{x,y)9, + P2{x,y)9y] Fi{A; x,y) , (76) 
Six,y)Fj{A + m;x,y) = [Qoix,y) + Qi{x,y)e^ + Q2ix,y)ey + Qs{x,y)ejy]Fj{A;x,y) , 

(77) 

— * 

where j = 2,3,4, m is a set of integers, A is a set of parameters, R,S,Pi,Qi are some 
polynomials and 9^ = xd^iOy = ydy). 

However, there is a special subset of values of parameters when the results of differen- 
tial reduction, Eqs. (1761) and ( 1771) . have a simpler form. This set of exceptional values of 
parameters can be defined from (i) condition that the hypergeometric function entering in 
the l.h.s. of Eqs. (il])-(ll3]), ([52])-(l5l]), ([63])-([65]), ([74]), ([75]), is expressible in terms of simpler 
hypergeometric functions (e.g. Gauss hypergeometric functions); (ii), condition that some 
of the coefficients entering in the inverse differential relations are equal to zero (infinity). 

For Appell's hypergeometric functions, Fi, F2, F3 and F4, the exceptional set of parame- 
ters are enumerated in Table [2J 

It is not surprising that the set of exceptional values of the parameters coincides with 
the set of parameters defining the condition of irreducibility of monodromy group of the 
corresponding hypergeometric functions (see [23] and references therein). We are expecting 
that the condition of irreducibility of the Mellin-Barnes integral pS] can be used as a criterion 
of irreducibility of Feynman Diagram. 

4 pfq - differential reduction of hypergeometric func- 
tion pFp—i. 

4.1 Non-exceptional values of parameters 

In this section, we will present the Mathematical^] based package pfq for differential reduction 
of the hypergeometric function pFp_i. In contrast to the version presented in [22] the current 

^It was tested for Mathematica 7.0. 
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version deals with non-exceptional and exceptional values of parameters. The Takayama 
algorithm is implemented in file "pfq.m", and an example of application is given in the file 
"example-pfq.m" . 

The program could be loaded in a standard way: 

<< "pfq.m" 

and includes two routines: ToGroebnerBasis[. . . ] and explicitForm[. . . ]. 

The main routine 

ToGroebnerBasis [Argumentsvector] , (78) 

calculates explicitly the ratio of the functions {Rk} and 5* of Eq. f fT9|) . The "Argumentsvec- 
tor" in Eq. flTSl) is the set of parameters of the hypergeometric function in the l.h.s of Eq. ([T9|) 
with an explicit set of integer numbers, by which each parameter should be shifted: 

Argumentsvector = {{a + M}, {b + K}, {x}} , (79) 

where M, K are integers, a, b are any symbols, and x denotes the argument of the hyperge- 
ometric function. The output of ToGroebnerBasis [. . . ] has the following structure: 

{{Qi, • • ■ , Qp+i}, {{ai, ■ ■ ■ , flp+i}, {1 + 6i, ■ ■ ■ , 1 + bp}, x}, factor}, (80) 

where 

• {{oi, ■ ■ • , ttp+i}, {1 + 6i, . . . , 1 -|- bp}, x} are the set of parameters and the argument of 
the resulting hypergeometric function entering in r.h.s. of Eq. f|T9l) : 

• {Qi: ■ ■ ■ ! Qp+i} are the rational functions where the numbering corresponds to power 
oi 9\ i = 0, . . . ,p; 

• factor is an overall factor. 

The explicit form of Eqs. (178l) -(l80l) is the following: 



/ ai + Mi, ■ ■ ■ ,ap+i + Mp+i 
b, + K,,--- ,bp+Kp 



X 



factor X (^Q, + Q,e + - ■ ■ + Qp+^9^y+^Fp 1+5/ .' ' .' ^) 



1. Example 1 : reduction of 3F2 

ToGroebnerBasis [ {{H-ai,2+a2, a3},{l+bi, b2+2},x} 

IntegerPart={l,2,0,l,2} change Vector={-1,-2,0,0,-1} , 



^All functions in the package HYPERDIRE generate output without additional simphfication. This 
is done for the maximum efficiency of algorithm. To get a simpler form of output, we recommended 
to use additional Simplify command. In particular, all considered examples correspond to FullSim- 
plify [ToGroebnerBasis [. . . ]]. 
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(a2-a3 + l)(fc2 + l) {b2+l)(xa'^-{a:iX-x+bi)a2+xai{a2-a3+l)+bib2) (b2+l)(-a3a;+x+(x-l)a2+fe2) 
(a2+l)(a3-fe2-l)' 



xaia2{a2+l){a3-b2-l) 



xaia2(a2 + l)(a3— 62-1) 



{{ai, 02, ag} , {61 + 1, 62 + 1} , a;} , 1} 



corresponds to 



3-^2 



ai + 1, 02 + 2, as 

h + l,b2 + 2 



(02 -03 + 1) (&2 + 1) 
(a2 + 1) (as - 62 - 1' 



(62 + 1) (a;a2 - {a^x - x + bi) a2 + xai (02 - as + 1) + 61&2) 



xaia2 (a2 + 1) (as - 62 
(62 + 1) {-asx + X + {x - 1)02 + ^2; 



xaia2 (a2 + 1) (as - 62 - 1) 



-1) 
3-^2 



ai, a2, as 
61 + 1,62 + 1 



2. Example 2: reduction of 4FS 

ToGroebnerBasis [ {{l+ai,l+a2, as,a4},{l+bi, b2+l,bs},x} 

IntegerPart={l, 1,0,0, 1,1,0} change Vector={-l,-l,0,0,0,0,l} 



ll — + - + — ^ 

' 02 63 ai ' ( 



f82) 



+02+^3 1 

aia2b3 ' aia2b 



■ j , {{ai,a2,as,a4} , {hi + 1,62 + 1,6s + '^} ,x} ,1} 



corresponds to 



4-rs 



1 + ai, 1 + a2, as, a4 
1 + 61, 1 + 62, 6s 

ai + a2 + 6s g2 _^ 1 ^ 
aia263 aia263 



4^3 



,1 1 1 
1+ I - + + - 

a2 6s ai 

ai, a2, as, a4 
6i + l,62 + l,6s + l 



X 



^3) 



In both these examples, IntegerPart = {■ ■ ■ } means the set of integer values of parameters 
in the original hypergeometric function, e.g. IntegerPart = {1,2,0, 1,2} in Example 1, and 
change Vector = {■■■} corresponds to values of the parameters which should be changed, 
e.g. change Vector = { — 1, —2, 0, 0, —1} in Example 2. 

Routine explicitForm[. . .] convert the results of reduction, Eq. ( IHTj) . to Mathematica- 
standard expression for generalized hypergeometric function. 

Example 3 : reduction of 2-^1- 

answer=ToGroebnerBasis [{l+ai,l+a2},{l+6i},x} ]; 



Integer Part= { 1,1,1} change Vector= {-1,-1,0} 
explicitForm[answer] 
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HypergeometricPFQ({ai ,a2},{6i + l},2;) , aia2a;(ai+a2— fei) HypergeometricPFQ({ai + l,a2+l},{6i+2},x) 
1—x {bi + l){aia2—aia2x) 

For non-exceptional values of the parameters the differential reduction is performed with 
the help of Eqs. (fT0|) . (fTT]) .( !T3|) .(!T ^ . The higher powers of the operators 9^ are expressed 
with the help of the differential equation for the hypergeometric function, Eq. ([9]). Also, the 
following relation is used in some cases: 



{Oj +m}p \ _ 11/0=1 IK^Kjrnj I ^ \ p ( L^-Jjp \ /o/fN 



m=l{K)n.} \dz) ' " \ {h}. 



4.2 Exceptional values of parameters 

When some of the upper parameters of the initial hypergeometric function are integer, then 
the higher power of the differential operator can be excluded with the help of Eq. fl2Tl) . 
In this case, the input of ToGroebnerBasis[. . . ] does not change, and the output of 
ToGroebnerBasis[. . .] has the following structure (case when only one upper parameter 
is integer): 

{{{Qi, ■ ■ ■ ,Qp}, {{l + oi, . . . , l + flp}, {2+6i, ■ ■ ■ ,2 + 6p},x},factorl}, {Qp+i, {}, factor2}}, 

(85) 

where 

• {{1 + ai, ■ ■ ■ , 1 + Op+i}, {2 + 6i, . . . , 2 + x} are the set of parameters and argument 
of the resulting hypergeometric function; 

• if upper parameter is integer, the proper is equal to zero; 

• {Qii ■ ■ ■ ) Qp} are the rational functions where the numbering corresponds to power of 
^\z = 0,...,p-l; 

• factor is overall factor; 

• Qp+i is the resulting polynomial entering in the r.h.s. of Eq. (1^ . 
In the explicit form 



/ Ml, ■ ■ ■ ,M^,ar+i + M^+i, • • ■ ,ap+i + Mp+i 



X 



factorl X ( Q, + Q2e + - ■ ■ + g/^-Mp+iF„ M' " " " ' ^ + ■ ■ ■ , l + a^+i 



2+6i,--- ,2 + 6p ^ 
+factor2 x Qp+i , (86) 

where M is a set of integers. 

Example 4 : reduction of 3F2 with an integer parameter. 
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ToGroebnerBasis [ {{3,l+a2,l+a3},{2+bi,2+b2},x} ]; 



IntegerPart={3, 1,1,2,2} change Vector={-2,0,0,0,0} 



(bi + l)(&2 + l)-x(a2+l)(a3+l) 
2{x-l) 



+ 1 



' 2 



1 / -x(a2+a3+2)+bi+b2+2 



x-l 



+ 3 



)} 



{{1, a2 + 1, as + 1} , {6i + 2, 62 + 2} , x} , 1 !>, { 

In an explicit form it is 

3, 02 + 1, cts + l 

2,^2 



(bi+l)(fc2+l) 



{},!} 



+ 



61 + 2,62 + 2 
1 / -X (02 + 03 + 2) + 61 + 62 + 2 



(61 + 1) (62 + l)-x (02 + 1) (aa + l) 
2{x - 1) 



2 V x-l 
(61 + 1) (62 + 1) 



3 ^ 



3-^2 



1, 02 + 1, «3 + 1 
61 + 2, 62 + 2 



2(x - 1) 

Other useful relations which are implemented in the function ToGroebnerBasis [. . 
are the relations derived in [26] : 



p-Tq 



61 + mi, ■ ■ ■ ,6„ + m„,a„+i 

61, ■ ■ ■ , 6n; 6„+i, ■ ■ ■ ,bq 



, • • • , dp 



ii=o i„=o 



where m,- are positive integers, Jn = ji + ■ ■ ■ + jn, and 



-4(ji,---j„) 



mi 
Ji 



^7) 



mn\ (62+W^2)ji(63 + "^3)j2 " " " (^n + J„-i (fln+l ) Jn " " " (Qp)jn 
(61) Ji (62)72 • ■ ■ {bn)jAb^l)jn ■ ■ ■ {K)jn 



Example 5 : reduction of 3F2 with integer difference values of parameters. 
ToGroebnerBasis [{{3+bi,l+a2,l+a3},{2+bi,2+b2},x}]; 

IntegerPart = {3 , 1,1,2,2} changeVector = {- 1 ,-1,-1} 

{( - ^ '^^t,\,^ , - ("2^+°3^-fel^-^+''l-fc2 + l)(fe2 + l) I JJq Q3} Jfe^ + 1} a;} 1} 
L 1 (x— l)(fei+2) ' {x-~ljxa2a3{bi+2) J5LL^) jjjl^' J5 J) J 
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In an explicit form 
3F2 



3 + 61, 1 + a2, 1 + «3 
2 + 61,2 + 62 



X 



^9) 



1 



(x-l)(6i+2) 
The following relations: 



(a2X + a3a;-6ix-a;+6i-62 + l) (62 + I) 
(x — I)xa2a3 (61 + 2) 



9 



02, ^3 
62 + 1 



X 



A, a 
l + A,b 



k rp 



{1 + A}„6 



A, a 

1 + 



i=o 



j=0 

{1 + A},_,,6 



p-L q 



1; 

2, {6fc},_i 



^n?:;(%-i) L 



{ttj — l}p-l 



(90) 



where g < r and a^, 6^ 7^ 1 are implemented in Mathematica 7.0. 



5 AppellFlF4 - Mathematica based program for dif- 
ferential reduction of Appell's functions ^1,^2,^3,^4 

In this section, we will present the mathematicaI based program AppellFlF4 for dif- 
ferential reduction of the Appell hypergeo metric functions Fi, F2, F3 and F4. The program is 
available in |38]. The current version, 1.0, deals with non-exceptional values of a parameters 
only. 

The program could be loaded in a standard way: 

« "AppellFlF4.m" 

The package includes the following basic routines: 

FllndexChange [changing Vector, parameter Vector] , (91) 

F2IndexChange [changing Vector, parameter Vector] , (92) 

FSIndexChange [changing Vector, parameter Vector] , (93) 

F4IndexChange [changing Vector, parameter Vector] . (94) 

and 

explicitFormFl[. . .] . (95) 
Hi was tested for MATHEMATICA 7.0. 
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The " changing Vector" in Eq. (!9Tl) - (|94l) is the set of integers at which we want to change 
the values of parameters of the Appell functions (vector fh in Eqs. (|76l) . ( |771) ). The set of 
parameters of Appell's function are defined in the list "parameter Vector" (corresponds to 
vector A + m and arguments x, y in the l.h.s. of Eqs. ( 176|) . ( 177|) ). We want to point out that 
enumeration of parameters and arguments in the " parameter Vector" correspond one-to-one 

to Eqs. dMD, dSD, dSSD, dSH). 

The output of FlIndexChange[] and F{234}IndexChange[], has a little different 
structure. 

5.1 Appell function Fi 

The structure of output of FlIndexChange[] is the following: 

{{A,B,C}, {parameter VectorNew}, {AppellFl}}, (96) 

where 

1. parameter VectorNew is the set of new parameters of Appell function Fi ; 

2. A, C are the rational functions corresponding to the ratio of Pq/R, Pi/R and P2/R 
of functions entering in Eq. ( 1761) . 

Example 6 0: reduction of Fi. 
FlIndexChange[{l,— 1,0,0}, {a,6i,62,c,^i,^2}] 

{ -fe^^feW' iStt} , + 1> - 1, h2, c, z„ Z2}, AppellFl} 

In an explicit form: 

Fi{a,bi,b2,c; Zi, Z2) = 
-azi + a + biZi + b2Z2 - c- zi + 1 (zi - l)(a - 61 + 1) Z2 - I 

7r w ~l ^2 

a — c + 1 (oi — l)(a — c+1) a — c+1 

xFi(a + l,6i-l,62,c;zi,Z2). (97) 

Routine explicitFormFl[- ■ ■ ] converts the result of reduction to Mathematica standard 
expression for the Appell function Fi. 

Example 7 : reduction of Fi. 

result = FlIndexChange[{l,— 1,0,0}, {a, 61,62, c,zi,Z2}] 
explicitFormFl [result] 

See footnote [2] 
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-^^^±i^^^(^^§^f^AppellFl {a + 2; h, h,- c + 1; z,, z,) 
^ (o+l^jgiziMAppellFl (a + 2; 61 - 1, 62 + 1; c + 1; Zi, z^) 

I (a(-zi)+a+b-izi+b2Z2-c-zi+l) 



a— c+1 



AppellFl (a + 1; 61 - 1, 62; c; zi, Z2) 



5.2 Appell functions F2, F3, F4 

The output of F2IndexChange[], F3IndexChange[], and F4IndexChange[], is similar 
and has the following structure: 

{{A,B,C,D}, {parameter VectorNew}, {NameOfFunction}, (98) 

where 

1. NameOfFunction is the name of Appell functions under reduction: AppellF2, AppellFS 
or AppellF4 ; 

2. parameter VectorNew is the set of new parameters of the Appell function; 

3. A,B,C,D are the rational functions corresponding to the ratio of Qo/ S,Qi/S,Q2/S 
and Q3/S of functions entering in Eq. f l77|) . 

Example 8 : reduction of F2. 
F2IndexChange[{0,0,l,l,0}, {a,bi,b2,ci,C2,zi,Z2}] 

1 Z2(a + zi(6l— ex)) 

az2{ci{zi-l)-bizi) I (zi-i)(62-c2+i) ci{zi-l){z2-l)-biziZ2 Z1+Z2-I 

Ci(zi-1)(62-C2 + 1) ' Cl ' Cl(2i-l)(fe2-C2 + l) ' Ci(^l-l)(b2-C2 + l) I ' 

{a, 61, &2 + 1, Cl + 1, C2, zi, Z2} , AppellF2} 
In an explicit form: 

F2{a,bi,b2,ci,C2,zi,Z2) = 

II ^\ U \ 1 Z2{a+zi(bi-ci)) 

aZ2 [Cl [Zi - 1) - biZi) ^ ^ ^ ^ (^i-l)(fe-C2 + l) g ^ Cl [Zi - 1) {Z2 - 1) - 0l^l^2 g 

Cl (2:1 - 1) (62 - C2 + 1) Cl Cl (2:1 - 1) (62 - C2 + 1) 



-2l + -22 — 1 

-6162 



Cl (Zl - 1) (62 - C2 + 1) 



^2(0, &i, &2 + 1, Cl + 1, C2; zi, Z2) ■ (99) 



The Appell functions F2, F3 and F4, are not built in the current version of Mathemat- 
ica (version 7.0), so that the AppellFlF4 package does not include a function similar to 
explicitFormFl [] . 
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6 Application: reduction of Feynman Diagrams 



6.1 Hypergeometric structure of Feynman Diagrams 

The standard way of analytical evaluation of multiloop Feynman Diagrams is based on ex- 
ploring of integration by parts relations [27] for the reduction of original diagram to set 
so-called master-integrals. This reduction can be performed in an arbitrary space-time di- 
mension (before construction of e-expansion) . Now there are a few ways of computer re- 
alization of this algorithm [28]. In a series of papers [l3l[Tl] it was pointed out that the 
hypergeometric representation of Feynman Diagram in conjunction with the differential re- 
duction algorithm allow one to make reduction directly in terms of hypergeometric functions. 
In particular, it was formulated the conjecture about the structure of the hypergeometric 
representation of Feynman Diagram (see as review [29]). 

Let us consider the hypergeometric representation of the Feynman diagram $(n, j; z), 

k 

$(n, J; z) = Y, Sain, J, z)H,0a, K f) , (100) 

a=l 

where Ha{/3a] Xa] the Horn-type hypergeometric functions defined by Eq. ([T]), j is a list 
of the powers of the propagators in the Feynman diagram, n is the space-time dimension, 
^ are the arguments of the hypergeometric functions which are related to the kinematic 
invariants of the Feynman diagram, {Pa,Xa} are the linear combinations of j and n with 
a polynomial coefficients, and Sa are the rational functions of the variables z with the 
coefficients depending on n and j. 

In general, the r.h.s. of Eq. f llOOp can include the different types of the hypergeometric 
functions (e.g. p+iFp and pFp_i, where p is integer, are allowed). Our proposition is that: 

(i) regardless the type of functions in the r.h.s. of Eq. fllOOp the number of basis elements 
is the same (up to a module of rational functions) for each term of the r.h.s. of Eq. f llOOP ; 

(ii) the number h of nontrivial master integrals following from integration-by-parts relations 
which are not expressible in terms of gamma functions is then equal to the number of basis 
elements L for each term of the r.h.s. of Eq. f llOOp . 

As illustration of the differential reduction algorithm and our conjecture, let us consider 
the diagrams shown in FigsUl El and |3l 

6.2 V()oo2 

Let us start from consideration of the two-loop propagator diagram V0002 depicted in Fig. ([1]) 
with one massive line and three-massless lines: 

rf"(fcifc2) 

[ik2-pyy'[ih - A;2)2]-i[A;2]-2[A;2_M2]i4 ' 

(101) 
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^0002 




■10002 



Figure 1: Diagrams implemented in package HYPERDIRE. Bold and thin lines correspond 
to massive and massless propagators, respectively. 



Such diagrams appear as one of the master-integrals in pole mass of heavy quarks. The 
Mellin-Barnes representation of Eq. f llOip reads: 

I '}^~"'^~ ^ ^ 

VoOOA/(jl,Cr2,(Ti, j4;P^) 



X 



dt 



^(-^)^(t+J4)^(^^-C^l-(T2-J4-^)^(jl + (Tl + (T2+J4-^^ + ^) 



ry r((Ti+fr2+j4-t+t)r(f-ji-ai-a2-j4-t) 

Closing the contour of integration in Eq. f ll02p on the right, we obtain 

[^i-^vrt]' r(f-ji)r(f-ori)r(f-a2)r(ai+a2-f 



(102) 



VboOM(jl,Cri,Cr2,j4;P 
^2 \ H 



(p2yi^i^^4-- r(ji)r(ai)r(a2)r(j4)r(r2-ai-a2) 

Y r(j4)r(n-ai-a2)r(ji+ai+a2-^) 

1-0-1-0-2) 

2 



MV r(ai + a2-f)r(f-j 



+ 



^ p I j4, n-0-1-0-2, f + 1-0-1-0-2 P_ 

^ ri + l-ji-ai-a2, f-ji-ai-a2 M2 

^■^■^^'"^''■^■"r(ji)r(n-ji-ai-a2)r(ji+ai + a2+j4-ri) 



M2 
X 3^2 



r(f)r(t-ji 

ji, ji + o-i + cr2+j4-n, 1+ji-f 
l+ji + c^i + cr2-^, ^ 




(103) 



Both hypergeometric functions in Eq. f llOSp are reducible to the 2-^1 functions with one unity 
upper parameter, namely. 



2Fl 



l,f+/i 



2F1 



/2 + ? 



(104) 



respectively, plus rational functions of z. In accordance with Section 16.11 there is only 
one nontrivial master-integral in the framework of integration-by-parts relations and a few 
integrals expressible in terms of the product of F-functions. 
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Figure 2: Diagrams implemented in the package HYPERDIRE. Bold and thin lines 
correspond to massive and massless propagators, respectively. 



In accordance with Tarasov's algorithm [30] there is one master integral of the same 
topology and two integrals expressible in terms of the product of F-functions [3T] . 

For completeness, we will present a few first terms of e-expansion of the master-integral. 
For n = 4 — 2e, where e is the parameter of dimensional regularization [32] and ji = ox = 
(Ti = (T2 = 1, we have 



Vooom(1, 1, 1, 1;"^ 



iFx 



1, l-£ 
3-3£ 



M2 



2l^-^ 



-2e 



-p 



M2 



F=^(l-£)F(l + 2£) p" 
(l-2e)F(3-3e) M2 

F(l+£)F(3-3£) 
Til-e)Ti2-e) 



1, e 
2-e 




^ (A05) 



where we use i ""{p 



2\n/2 



-p 



2\n/2 



. For the second 2-^1 hypergeometric function the all-order 
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Figure 3: Diagrams implemented in the the package HYPERDIRE. Bold and thin hnes 
correspond to massive and massless propagators, respectively. 



£-expansion are available: 

1-s 



2F1 



1, e 
2-e 



2{l-2e)z 



i=i k=i 



(106) 

where Li„ (z) and ^(z) are classical and Nielsen polylogartihms, respectively. For the first 
hypergeometric function we apply the following relation: 



(1 - 2g) 
2-3£ 



1, 1 - £ 

3-3£ 



;i - z) 



1, l-£ 

2-3£ 



(107) 
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The first coefficients of tfie e expansion read: 

r_ _ 2 .2 ^)00^/(1.1.1,1;^') ^ 1 T\l-e)T{l + 2e) 

+£<(! + ln(l - 2) + 2 In(-z) 



1-z 



+e^3 + 2\n\-z) +4(2-2 [\n\l-z)-\n{-z) ln(l-2) + 2Li2 (z)] 



+8^9 + ^ ln\-z) + 8C3 + 8C2 ln(-2) 



+2- 



1 - z 



4Si,2 (2) + - ln3(l +4 ln(l -2)Li2 (2) -SLig (z) 



+2C2ln(l-2)-ln(-2)Li2 (z)-ln(-z) ln2(l-2)+ln^(-z) ln(l-2) 



+e^{ 27 + ^ ln\-z) + 40(4 + 8C2ln'(-2) + 16C3ln(-z) 



1 - z 



881,3(2) - 1082,2(2) + 14Li4 (z) + \ \n\l -z)+ 2C2Li2 (z) 



+8 ln(l - 2)81,2 (z) - 10 ln(l - z)U3 (z) - 4(3 ln(l - z) 
+41n2(l -2)Li2(2) + 2C2ln^(l -2) +ln^(-2)ln^(l -2) 

-2 ln(-2)8i,2 (z) + ln(-2)Li3 (z) - \ \xv{-z) - z) 
-4C2 In(-z) ln(l -z)-2 \^{-z) ln(l - z)Li2 (2) 



+ ln2(-2)Li2 (2) - - ln=^(-z) ln(l - z) 
o 



(108) 



wliere z = p^/M'^. 

To cfieck Eq. (I108p . we evaluate tlie ffist few coefficients of tlie e expansion of tlie original 
diagram in the large-mass limit [33] using the program packages developed in Refs. |34] to 
find agreement. 



6.3 £"10002 

Let us consider the 3-loop bubble diagram -E10002 depicted in Fig. ([T]) with two massive lines 
and three massless ones: 



^iooo2(A, ai, a2, (3, a, m^, M^) = j — 



- ^3 ) 2] «i [A;2] «2 [( - ) 2] /3 [A;2 _ ^2] A [ A;2 - M2] - ■ 

(109) 
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Its Mellin-Barnes representation reads: 

r(f-«i)r(f-a2)r(f-/3)r(ai,2-t) 



X 



X 



r(ai)r(a2)r((T)r(A)r(/3)r (§) v{n-a^; 



dt 



m' 



y r(-t)r(/3+t)r (a;3,,-f +t) r (n-ai,2+t) r (ai,2,A-?^-t) 



M2 / 



r(t+f) 



(110) 



where 



«i,2,/3,A,a = ai + a2 + \ + (3 + a , 01,2, a,/3 = ai + a2 + X + (3 . 
Closing the contour of integration in Eq. f lllOp on the right, we obtain 



(111) 



X 



^looo2(A,al,a2,/3,f^,"^^M2) = [i^ "tt^]' ^ 

r(f-^i)r(f-a2)r(f-/3)r(a,,2-t) 
r(ai)r(a2)r(a)r(A)r(/3)r (f ) r(n-ai,2) 

r(/3)r(a^,;3-t)r(n-ai,2)r(ai,2,A- 

r(f) 

m2V^'^'^-"r(A)r(ai,2,A,/3-n)r(ai 



-m 



2-\n-ai,2,A ('_]\^2-\n/2-a^,, 



3-^2 



Z?, ci/3,cr — ^, n — ai 2 



f,l+n 



+ 



'^l,2,A 

2,is,\,a-^) r(n-ai,2,A) 



M2y 



M2 y 

X3i^2 



r(ai,2,A-t) 

A, ai,2,A,^ — ai,2,A,/3,o- — 
ai,2,A-t,l + ai,2,A M2y 



(112) 



Both hypergeometric functions in Eq. f lll2p are reducible to the 2-F1 functions with one unity 
upper parameter, namely. 



n 
2 



"+/2 



1,/r 



3n 
/2-lf 



(113) 



respectively, plus rational functions of z. For this diagram there is only one nontrivial 
master-integral in the framework of integration-by-parts relations and an integral expressible 
in terms of the product of F-functions. 

For completeness, we will present a few first terms of the e-expansion of the master- 
integral. For n = A — 26, where e is the parameter of dimensional regularization [32] and 
ai = a2 = cr = A = /3 = l, we have 

£;looo2(l,l,l,l,l,m^M2) = [z7r2--]'^-2M-2s 



[m 



F(l-£)F2(l + e)F(l + 2£) 

2£3(1 -£)2(1 _2£) 



2F1 



1,6 
2-6 



-l + 2£ 



F(2-£)F(l + 3e) 
3(l-3£)F2(l + e) 



( l,-l+3£ 


rn? \ 







(114) 
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The all-order £-expansion for the first hypergeometric function is given by Eq. (11061) and for 
the second hypergeometric function we apply the following relation: 



2F1 



l,-l + 3£ 

l+e 



l-3£ 3£(l-3£)^(l-z) / l,l + 3£ 
l-2e^^ {l+e){l-2e) ' ' V 2+£ 



The first coefficients of the e expansion read: 
Eiooo2(l,l,l,l,l,m2,M2^ 



T{1- e)T{l + 2e)z 

(l-£)(l-2£)2r(l + £) 



X 



+ 



^(l + 2z) + ^[l + 2(l-z)ln^] 



3e 



In z + (1 - 4z)C2 + (1 - 2) (In^ z - 3 In z ln(l - z) - 3Li2 (z)) 



3 1 1 fin 2 2 

+ 2 + 3C2 + 3C3 + ln^ + -ln2^-(l-^)|yC3-gln'^-2Si,2(^)--C2ln^ 

+ ln(l - z) [2C2 - 2Li2 {z) + In^ 2] - In z In^ (1 - + In ZU2 {z) | 
+^ <! ^ + C2 - ^Cs - + ^ In^ 2; + In^ z + i (9 + 2(2) In ^ 



+ {l-z){ 7C4-Li4 {z)+\nz 
1 



2Si,2(^)+Li3 W--C3 



+^ln^^-4Si,3(;2) + ln(l-z) 

y 



4C3-4Si,2(^)--ln3z 



+21n(l - z) 



ln(l — z) — In z 



C2 - Li2 {z) 



+ ln^;zln^(l-2) 



— lnzln^(l - z) -ln^2; 
3 ^ ' 









Li2 {z) - 




111 



(115) 



(116) 



where z = rn? /M'^. For z = 1, Eq. (lll6p coincides with the results of [35]. 
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6-4 -B?i22000 ^ -^112200 ^ -^11220 ^ ^1122 

Let us consider the g-loop bubble type diagram -B1122000 E (see Fig. E]) defined as 

^112200o("^^ "1' "2, /32, CTl, ■ ■ ■ , (^x, Pi, ■ ■ ■ , Pr, 1^1, ■■ " , l^z) 

■ ■ ■ fcg) 



Pr 



X 



21<7i 



where P = kr+x+i + ■ ■ ■ + kg . and r, s and x are the number of massless fines. Tfie Mefiin- 
Barnes representation for diagram -B1122000 

^1122000("^^ «b «2, /32, CTl, ■ • • , CT^, Pi, • • • , Pr, 1^1, " " " , l^s) 



^_^2^ 2-^1,2 (^_^2^ 2 W"^>^P>-'''>ft,2 



r(ai)r(«2)r(/3i)r(/32)r (f ) 



X 



X 



r(cr.) J I Tip,) i I TK) i r(fr-p) 
/ rftf r(a-fx-t)r(a;3i,fe,.,p,a-t(g-i)-t) 

r(ai+t)r(a2+t)r(ai,2-t+t)r(f+t)r(f(r+a;)-ap,<,+t) 



r(ai,2+2t)r(f(x+i)-(T+t) 



r(a^,,-f(r+a;-l)-t)r(/3i,2 + 2a,,^,,-n(g-2)-2t) ' ^ ^ 

where x + r + s = q — 2,q>3, and we have introduced the short notation, 

r X s 

a=l a=l a=l 

ai,2 = ai+a2, /3i,2 = /3i + /32 , Oi/.p.a = cr+'^+P , ctp.fT = cr + P , 

a/3i,i^,p,(7 = /3i + o- + '^+P , a/32,i^,p,a = /32 + o-+«^+P • (119) 

The resufi could be written as a sum of 5 hypergeometric functions of type gFy of argument 
m^/M^ or M"^ /m?. The expression is really long to be considered here. 

More interesting particular cases could be derived when g, r, s, x take special values. 

Diagram m Fig. a emerges from with . = 0, [jllj-^] = 1 and 

x + r = g— 2. In this case, = and Eq. (11181) is simplified due to limi,_j.o — , "'''''^ n^, — -— = 1 . 

r(ap,,j— 2 ((?-3)-s) 



^We are indebted to A. Grozin who suggest us apply our technology to analyze diagram of this type. 
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The detailed analysis of this diagram was presented in [I3] . The result for diagram -B112200 '^^^ 
be expressed as a sum of four hypergeometric functions of type yFg or of lower order function. 
The general result again is to be published. We present only the results of differential 
reduction. In the case when the powers of the propagators are all integers, the appearing 
hypergeometric functions are reducible to the following hypergeometric functions and their 
derivatives]^ 



{1,9,9^,9^} X 5F4 



l),/2-fx,/3-f(x + l),i + /4 + f(g-x-2) 



/5 + f,/6 + f(g-x-l),i+/7-fx 

i + Ji-f,/2-f(g-2),/3-f(g-l),J4-fg 
/5-t(g-2;-l),/6-t(g-a;-2),i+/7-f(g-l 

l,i+/i,/2-f(g-l),/3-t(g-2),/4-f(g-3) 
/5 + t,|+/6-t(g-2),/7-t(g-x-2),J8-f(g-x-3) 



(120) 



where Ik are arbitrary integers. In the last expression, some rational functions are also 
generated. In accordance with Section lGTTI this integral has four master-integrals for arbitrary 
values of x and q. Further simplification is possible for (i) x = and (ii) x = 1. In these 
cases we have: 



for X = 



{1,^}X 

{1,9} X 2F1 
{l,^}x 3F2 



1,/. 



2-f,|+/3 + f(g-2) 
/4 + f(g-l),i+/5 



1 _ n h-Hn 

2 2 ' ^ 2^ 

^3 + i-t(g-i) 



l,|+/i,/2-f(g-l) 
/3+f./4+|-t(g-2) 



(121) 



and in accordance with Section 16. H this integral has two master- integrals only, 
for X = 1 



{1,9,9^} X 4F3 
{1,9,9^} X 3F2 
{1,9,9^} X 4F3 



l,/i-f,/2-n,i+/3 + f(g-3) 
/4+f,/5 + f(g-2),|+/6-t 

/i+i-f,/2-t(g-i),/3-tg 
/4-f(g-3),/5 + i + f(g-l) 

l,i+/i,/2-f(g-l),/3-f(g-2) 
/4+i,/5+|-i(g-2),/6+i(g-4) 



(122) 



and in accordance with Section 16. H this integral has three nontrivial master-integrals only. 

For a next diagram, -B11220 (^^^ -^^S- HI)' the Mellin-Barnes representation follows from 
Eq. f lllSp putting r = 1, p = pi — )• 0, 2; = and x = q — 3. For simplicity, we will present it 



^Two of the original hypergeometric functions have a similar parameter structure and the same basis 
functions. 
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an explicit form: 



^1122o("^^ "1' "2, /32, CTl, 

n "1 <7 ?i n , , , 



n9-3^(2 



X 



r(aOr(«2)r(/3i)r(/32)r (f ) r"=^ r(a, 

1 rft(^) r(a-|(g-3)-t)r(a;3,^,.-f(g-l)-t)r(f+t) 



r(ai,2 + 2t)r (/3i,2+2a-n(g-2)-2t) 
and in agreement with Eq. f lllQp we have introduced the short notation, 



(123) 



(124) 



The result could be written as a sum of four hypergeometric functions of type qF^ of the ar- 
gument m? /M'^ or M"^ /iv?. Again result is to long to be published. The results of differential 
reduction follows from Eq. (I120p by substitution x = g — 3: 



/i-|(g-4),/2-f(g-3),/3-f(g-2),l + /4 + f 
J5 + f,/6+n,i + /7-f(g-3) 

i+/i-|,/2-f(g-2),/3-f(g-l),/4-fg 
/5-n,/6-f,i + /r-f(g-l) 

i+/i,/2-f(g-l),/3-f(g-2),/4-f(g-3) 

J5-|,/6+f,K^7-i(g-2) 



(125) 



Further simplification is possible for (i) g = 3 and (ii) g = 4. In these cases, we have: 
• for g = 3 , 



{1,^^}X 3i^2 
{1,^}X 2F1 
{1,0} X 3F2 



-'-,-'2-25 2+-'3+2 

h+n, I + /5 

-^1 + i ~ f ,-^2-|^ 

/3 + l-n 
1, | + /i,/2-n 

r I n Til " 
i3+2'-'4+2~ 2 



(126) 



and in accordance with Section 16. H this integral has two master- integrals only 
for g = 4 



{i,e,e^}x 4F3 
{i,e,e^}x 3F2 
{i,e,e^}x 3F2 



l,/l-f,/2-n,i + /3 + f 
/4 + ^,/5 + ?^,H/6- 



14-1 2 ' "0 ' 2 
-^1 + I — f 5 -^2 — 1?^, -^3 



2'-, -'3 
2 

i+/i,/2-|n,/3-n 



2 
2n 



(127) 
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and in accordance with Section 16.11 this integral has three nontrivial master-integrals only 
(see [36]). 

The three-loop diagram Bf^22 depicted in Fig. corresponds to a particular case of 
diagram -B11220 for a = and q = 3. In this case, n^~3 r(n/2-g-o _ g closing the contour 
of integration on the left we obtain: 



Sl\22("^^^^al,a2,/3l,/32) 



n - 


3 n 


zi-"7r2 


(_^2) 2-01,2 ^_^2)n-/3i,2 







r(«i)r(a2)r(/3i)r(/32)r (f ) 
r (A-f) r (/32-t) r («i,2-f) r (f) r(ai)r(«2) 



r(ai,2) 

a, ao a, o-^ ^ t^L2±l t^l^ 

p I "Ij'-^ij'-'l.i 2' 2' 2 ' 2 

' + 1-/32 + 1, l-/3i,2+n,^,^ 



r (f-/gi) r (f-/32) r (/3i,2-t) r {aa„p,,,-n) r (aa,,;^,,,-^^) r (a,^^^^^^-f ) 

r(ai,2+2/3i,2-2n) 

X ( — — I F I "l.2~2'^"i./3i,2~'^'"a2,/3i,2~'^5'^ai,2,/3i,2~T' 2 ' 2 

V^V ' ' V l + /3i-t,l + /32-t,l + /3i,2-n,^^ + /3i,2-n,^^if±i + /3i,2-n 

r(/3i)r (/32-t) r (f r r (a2 + /3i-t) r («i,2 + /3i-n) 



2\ /3i-n/2 



6-^5 



r(ai,2 + 2/3i-n) 
/3i,ai-h/3i-f ,a2 + /3l-?,al,2+A-?^, 14 



/3i-/32 l+/3i-/32 



2' 



l+/3i-/32,l+/3i-^, 1-/32 + 



2 ' 2 



r(/32)r (/3i-f ) r {^-/s^) r («i+/32-t) r (a2+/32-t) r K2+/32- 



2N /32-n/2 



r(ai,2 + 2/32-n) 
/32,ai + /32-f,«2 + /32-t,«i,2 + /32-n,^±%^ 



l+/32-/3i, 1+/32- f, l-/3i + t, ^+/32, 2 



ai,2+l-^ I 



(32 



(128) 



6.5 



J'i kr J'i 



Let us consider the g-loop propagator type diagram Jq23 (see Fig. [3]) defined as 
^023(^1 ^^2^«b«2,c^l, 



[kl - M2]-i [fci - M|]"^ [fci]- ■ ■ ■ [(p - A;i 



O'q-l 



(129) 
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The Mellin-Barnes representation of Eq. f ll29p is 




dtidt2T{-ti)T{~t 



r(ai)r(a2 



'-^1=1 



Ml 

p2 



2 \ *2 



r (|-ai-ti) r (1-02-^2) r {aa,,a^,a-n+ti+t2) 

r (|^^ — C[cti,02,o-~^1~^2) 



(130) 



where 



q-l 



a 



cr,; , a. 



2 , ^CKi,CK2,cr 



a2+cr 



(131) 



j=i 



Closing the contour of integration in Eq. fll30p on the left, we obtain 



0^023(^1 5 ^2^ "2, C^i, ■ ■ ■ , Cr^-i) 



r(ai)r(a2)r (f ) 



ir(t-^0 



X F4 I a, 



2' 772 



2\ ^-"2 



+ (4) ^ ("'^ - 2) ^ 2'« - ^)) ^ 2'«-2) 

/ft Tl Th 

X F4 (^flai.a - -(g - 1), - 2 (^~2), -, l-a2 + 



f_ Ml 



(132) 



Let us analyze the result of the differential reduction of hypergeometric functions in Eq. fll32p 
assuming that all powers of propagators are integer. For q = 2, both hypergeometric func- 
tions are reducible to function and their first derivative only (in accordance with the results 
of Section 13. 6p 



F4{a + mi, 6 + Ji, c + m2, b\x, y) 
Fi{a + mi, Ji, c + m2, h + msjx, y) 



:i, dy) X F4(a, 6, c, 6|a;, y) + i?i(a;, y) , (133) 
:i, 9^, 9y) X F4(a, 1, c, y) + Riix, y) , (134) 



where {mj} is a set of integers, and Ri, R2 are rational functions. For g > 3 the result of 
differential reduction coincides with Eq. (I77|) . 

In accordance with Section 16. H there are three nontrivial master integrals for q = 2 and 
four ones for q > 3. 
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For completeness, we also present the explicit results for the diagram Jq^2 (ses Fig. [3]) ( 
when external momenta on mass shell of one of the internal masses) 



JJl2("^^ a, A, ai, ■ ■ ■ , ag-i] 



[kl - m?Y[kl - M^Y[klY^ ■■■[{j>-ki- ■■ 
The hypergeometric representation for this diagram is 

n 



k.fY"-- 



(135) 



r(a)r(A)r (f ) 



_ir(f-a,:) 



r(cr.) 



X 4F3 



'-*'a,A,(T '^i'^ 2^^ 2' 2 



X I 2 7 ^ ' 2 ' 



4m^ 
M2" 



n 



X 4F3 



a,,,-f(g-l),a-f(g-2),^ 



n— g n— a+l 
2 ' 2 



4m^ 



(136) 



where 



9-1 



(137) 



In the case when the powers of propagators are all integer, the result of the differential 
reduction of the hypergeometric functions in Eq. (1136^ are the following (see also discussion 
in [37]): 



for q = 2 



{1,9} X 2F1 
{1,0} X 3F2 



t + /3 



In accordance with Eq. fl6.ip . there are two nontrivial master-integrals only. 

{i,e,e^}x 3F2 
{i,e,e^}x 3F2 



:i38) 



for g>3 ^2. . f h-lq,\ + h,h-l{q-l) 

/i-f(g-2),/2-f(g-l),i + f + /3 



2 2 ' 2 

n + /4, f + Is 



(139) 



In accordance with Eq.f l6.ip . there are three nontrivial master- integrals only. 
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7 Conclusion 



The differential reduction algorithm [9] allows one to compare Horn-type hypergeometric 
functions with parameters whose values differ by integer. In this paper, we have presented 
the Mathematica based package HYPERDIRE for differential reduction of the generalized 
hypergeometric functions p+i-Fp and Appell functions Fi, F2, F3 and F4 to a set of basis func- 
tions, defined by Eqs. (|T9l) . (1761) . ( 1771) . respectively. These functions are closely related with a 
large class of Feynman Diagrams [5] . In contrast to existing packages, our package performs 
reduction of hypergeometric functions before e-expansion and works with arbitrary values 
of parameters. The package could be easy extended to other Horn-type hypergeometric 
functions by adding new modules. 

As an illustration of our approach, we considered a few examples, namely the Feynman 
diagrams shown in Figs. [T]|3] with arbitrary powers of propagators and space-time dimension. 
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